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Abstract 

The negative dimensional integration method (NDIM) is a technique where several difficulties 
concerning loop integration can be overcome. From usual covariant gauges to complicated Coulomb 
gauge integrals, and even the trickiest light-cone integrals one can apply the methodology of NDIM. 
In this work we show how to construct a general formula — we mean arbitrary exponents of 
propagators, off-shell external momenta and distinct massive propagators — for one-loop scalar 
integrals, for covariant gauges, and apply it to one through six-point loop integrals. We present 
detailed analysis of pentagon and hexagon scalar integrals for massive/massless internal particles 
being external momenta on or off mass shell. 
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I. INTRODUCTION 



Important mathematical methods have been required to evaluation of the complex Feyn- 
man integrals in the calculations of scattering amplitudes in QED and QCD, in the radiative 
corrections, study of Green function behavior, renormalization group and others problems 
in quantum field theory. The integration using Mellin-Barnes representation |2], EL the 
Gegenbauer polynomial technique [|J, integration by parts ||, negative dimensional integra- 
tion (NDIM) |], [<J, string inspired methods||, differential equation approach|J and several 
others |I| |TT], 0, [IS], |I§ |TJ, |T7|, 0, 0, |2D|, [2TJ, are some of the technique that have 
been currently used. 

In the present article we make use of the NDIM || [/]. Its implementation transfer the 
complexities of the performing D-dimensional Feynman integrals to resolution of a system 
of linear algebraic equations that we call of system of constraint equations. Therefore, when 
we choose one solution to those constraint equations we are obtaining one specific solution 
to Feynman integral in question that is defined in a specific kinematic region. The others 
solutions to system constraint equations represent the possibilities of analytic continuation 
of such solution. In other words, the NDIM solve the Feynman integral as well as obtain 
solving the constraint system equation, all the analytic continuation possibilities to the 
solution found. We note that after arrive the final formula to Feynman integral is necessary 
yet carry out other but straightforward analytic continuation, now to negative value of the 
powers of the propagators in the scalar integral. 

We implement the NDIM in this paper to obtain the one-loop general massive n-point 
function. As already mentioned above we obtain, by solving the system constraints equation, 
all the other solutions analytically continued and show that the number of this possibilities is 
a function of number of external lines n. We present then applications to cases of 1 to 6-point 
functions, considering only scalar integrals since any tensorial integral can be reduced to 
scalar integral according to [22||. The hypergeometric series we choose to present in this work 



are such that can be used in the dimensional or analytic regularization schemes |23[ |2J, |25 



always preserving gauge symmetry, a fact well-known from quantum field theory. Similar 
results for N-point scalar integrals were obtained in |], |2|] using Mellin-Barnes approach, 
however Davydychev quote only the hypergeometric series he call "symmetric". Here we 
write such results and several others, an interesting feature of NDIM, providing a very large 
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number of hypergeometric series that represent the original Feynman integral. 

This work is organized as follows. In the early section we present a detailed approach to 
the implementation of the NDIM to one-loop n-point function. In the subsequent sections 
we have applied the results from section one, the general formula, starting from one-point 
integrals, two, three, four (box), pentagon and finally hexagon scalar integrals. These results 
are exact, i.e., no approximations was made. The solutions obtained in this sections are given 
in terms of the hypergeometric functions (see appendix) and compared with known ones, 
when they are available, in the literature. 

II. ONE-LOOP N-POINT FUNCTION: GENERAL FORMULA 

In this section we present the calculations to evaluate the one-loop scalar integrals with 
the NDIM. Consider a one-loop Feynman diagram with n, n = 1,2,..., internal momenta 
/ , Iq — h,lo — h, ---Jo — l n -i an d masses m , mi, m 2 , m n _i, where l±, l 2 , l n -i are given 
in terms of a linear combination of the external momenta. Its scalar integral associated is 



J i=i 
where cuo, on are positive parameters. Then, the exponential function above can be expanded 
and we get 




(1) 



Consider now the gaussian integral 




n-l 





x 



J d%[ll - mlp[(l - h) 2 - mlr...[(l - Z n _0 2 




(3) 



Using the definition 



n-l 



a = a + ^2 



(4) 



i=i 



we can rewrite the I integral of form 




a 



(5) 



i=i 



i=0 
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After evaluate the integral we have 



fyn-l Q,.7.y n-1 n-l 

" ^ D/2 exp{ ^ =1 ^ - £ ^ + £ «*"»?}. (6) 



a a i=l i=0 



Yet can us rewritten this expression of form 



a a 

n-l 



+ E«^I- (7) 



i=0 

Also can be show that 

n— 1 n— 1 re— 1 

E 'lj = 2 ^ ttittj-ii • + E ( 8 ) 

i,j = l i>j = l i=l 

n—l n—l n— 1 

E "i«i^ 2 = E + ?) + E ^ 2/ i • ( 9 ) 

ij'=l i>j=l i=l 

Performing the substitution of this results in (|35|), we get 

^re— 1 _ „ /2 _ r2 



, _ r - ... n; , i>j=l a i a 3 l ij ~ a 0z2i=l a i l i 

a, a 

n-l 

( 10 ) 

where = /j — Zj. From exponential argument above we have w = " 2 ~3"+ 2 terms with 
different coefficients namely 

OL\OLi-, «1«3, ai«4, ... 

OLlCLz, 02^4, ••• a2«n-l, 

0*30*4, ... a 3 a„-i, 



••• «n-20„-l. 

Also, there are n — l terms with different coefficients a^on and n terms to coefficients a». 
This result in w + 2n — 1 = n!±n terms with different coefficients. After expansion of the 
exponential above We can write 

oo -i 

J = (K )D /2 1 



Q/ . ^— ' Q>jl+32 + ---+jw + n-l 

jl,—Jw+2n-l 
v Jl+32+-+jn-lj'l +in +•»+ j2n-2 
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v ^h+in+j2n~2 + ---+hn-4, 32n-2+j3n-4,+j4n-T+--+jw 

xa 2 ...a n _i 

jn- jn+1- jw+n— 1- 

^ (m§) ,u, +" (mf) Jt0+n+1 (m^_ 1 )- 7u,+2n - 1 ^ 

jw+n- Jui+n+l' jw+2n— 1- 

If we take the multinomial expansion in the exponents of the a, given by (0), we get 
1 1 



Q,D/2+il+i2+-+iw+™-l "~" [q, _)_ ^^J 1 Q ; i ]-D/2+Jl+i2+-..+i«,+n-l 

£ T(l-D/2-j 1 -j 2 -...-j w+n ^) 



Jw+2nviJju+3n-l 



x a*__ai -, (12) 

3w+2n- 3w+2n+\- 3w+3n-l- 



with the constraint 

D/2 = —j\ — ]2 ~ ■■■ ~ jw+n-l ~~ jw+2n ~ jw+2n+l ~ ■■■ ~ jw+3n-l- (13) 



Performing the substitution of this expression in ( IT ) and compare the exponents of the 
parameters in (|3|), we have follow constraint equations 

OO = jl+j2 + -+ jn-1 + jw+n + jw+2n, (14) 

a l = jl + jn + ••• + J2n— 2 + Jw+n+l + Ju>+2n+l; (15) 

a 2 = J2 + jn + j2n-2 + ••• + j3n-6 + jiu+n+2 + jw+2n+2, (16) 

a 3 = jz + jn+1 + j2n-2 + jzn-5 + ••• + jln-10 + jw+n+3 + Jw+2n+3 (17) 



«n-l — jn-1 + j2n-2 + j3n-6 + J4n-10 + ••• 

••• + Jw+n— I ~\~ 3w+2n-l + jw+3n-l (18) 
D/2 = —ji — J2 — ... — Jto+n-l — jw+2n — jw+2n+l ~ ■■■ ~ jw+3n-l- (19) 

Using the results we have S = w + 3n — 1 = " 2 + 3n sums with n + 1 constraint equations 
and Cs, n +i different forms to evaluate. But, if there are Fj internal lines or F e external lines 
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associated to no massive fields, the number of the sums reduce to S — Fi — F e . Then, the 
number of the different forms to perform the S — Fi — F e sums is 

f n 2 +3n — p. — p\\ 

C S - Fi - Fe>n+1 = (n + 1)! 2 ( H!^2 l _ F 7T^j- ( 2 °) 

Finally we obtain the final expression to integral given in (D, that is 

j(n) = j( n )(D,a , ...,a n - 1 J 1 ,l 2 , ■■■,ln-i,m ,m 1 , ...rrin-i) 

= J d D k[ll - mlri(l - hf - m?r...[(Z - ln-if - rr&.J"- 1 

= 7r D / 2 (-i)£r>r(i + ao)r(l + ai)...r(l + a n _ x ) 
- T{l-D/2-j 1 -j 2 -...-j w+n „ 1 ) 



x 



E 



n ,...,j s Jw+2n ] -Jw+2n+V---JS- 



Jn-l 



3l ] J2i jn-r- 

X ; ; . . . , 



Jn- Jn+l- Jw+n—l- 

x {ml) jw+n (mf )^+"+! (m 2 _ 1 )^+ 2 - 1 

Jui+n' Jtu+n+1- 3w J r2n— 1- 

This expression only represent the one-loop n-point function after the analytic continuation 
in the parameters a , ai, a n _i to negative value. 

III. ONE-POINT FUNCTION 

One-point functions at one-loop level are the simplest Feynman loop integrals and we 
start with them for completeness. The integral associated to one-loop one-point function, 
case n — 1, given by 

J (1) {D, at, m)= J d D l (l 2 - m 2 ) a \ (22) 

can be evaluated by method above, obtain, after analytic continuation to i < 0, the known 
result 

J«(Aai,m) = / /2 (-«i)- B M-™ 2 ) ai+D/2 (23) 



that is to according to [26] when a\ = — 1, and the Pochhammer symbol is defined as, 



(a\b) = (a) b = ^±^, (24) 



we will turn to the left form when we deal with pentagon and hexagon integrals, because the 
number of sum indices will be large (more than 10 sometimes) and then the first notation 
becomes better to read. 



IV. TWO-POINT FUNCTION 

Two-point integrals are needed in order to study radiative corrections such as self-energy 
and vacuum polarization. These integrals raises no difficulties and their results are well- 
known of quantum field theory courses. The integral associated to one-loop two-point func- 
tion, case n = 2 given by 

J (2 \D,a 1 ,a 2 ,h,m ,m 1 ) = J d D l (l 2 - m 2 ) ai [(l - k) 2 - m 2 ]° 2 , (25) 

can be evaluated by method above, one obtains ten different solutions analytically continued 
that can be calculated by general expression 

, ai,a!2,Zi,mo,mi) 

= 7r D / 2 (-l) ai+Q2 r(l + a 1 )r(l + a 2 ) 

- T(l-D/2- 3l )(-l 2 ) K) (mf) 

h h=o (%(% (Ik (!)*(!)* ' 1 } 
using the constraints equations 

D/2 = -h-u-js, (27) 

on = J1+J2+J3, (28) 

Oil = JI+J3+J5- (29) 

The two point function will be obtained after the analytic continuation of the each solution 

to Cti, OL2 < 0. 

We choose one convenient solution given by (consider a 2 = a± + a 2 + D/2) 



,ai,« 2 ,li,m ,mi) 
= 7r D / 2 (-m2)^{(-« 2 )^ 1 -D/2p/2), 
xF 4 [-cr 2 , -«i; D/2, 1 - «i - D/2 

7 



+( 



,m 



2 . 

mf 



0\ ai +D/2 



■D/2 



/2 777^ 

xF 4 [-« 2 , D/2; D/2, 1 + a, + D/2\^; -§]}, 

m{ m{ 



(30) 

where _F 4 is an Appel's hypergeometric function ||27|| . This solution contains special cases. 
For example when mo = 0, mi = m this expression can represent the integral associated to 
self energy of fermions. In this case the solution above can be written of the form 



J(2) = j(2)p jaijQ , 2> / lj0j 



m) 



7T 



D/2 (-« 2 )-a 1 -D/ 2 P/2) Ql (-m 2 ) 

XaFif-^-aisD^I^], 



(31) 



where 2^1 is the Gauss hypergeometric function [T?\. The analytic continuation of this ex- 
pression can be obtained by choosing of others three indices between ji, j$. This solution 
is divergent only to great value of l\ and the wait infrared divergence not appear here. If 
we perform dimensional regularization with D = 4 — 2e, ol\ = a 2 = — 1 and e —>■ in the 
equation above, we get 



j(2) = j(2)(4_2e,-l,-l,Zi,0,m) 



A - 7r 2 [log(-m 2 
J 2 



rn 



I 2 



tog(l-^) + - r log(l--U-l]. 



m z 



77T 



(32) 



where A 



7T 



2/1 



7 — log7r). This result is to according to [2£|. Other case of interest 



we have when m-o = m\ = m. This case is associated with the diagram of the vacuum 
polarization. The expression (|30|) is then given by 



, a>i, 0:2, h, m, m) 



-02, — Oil, — «2 
2 ' 2 



I 2 



Am 2 



(33) 



where 3F2 is a hypergeometric function. The above solutions as well as its analytic contin- 
uation is to according to U and . 



S 



V. THREE-POINT FUNCTION 



We turn to scalar three-point functions, increasing a little bit the difficulties, we one put 
the external momenta off-shell and consider three propagators with distinct masses. The 
one-loop three-point function is given by 



, an, "2, oi 3 , h, l 2 , m , mi, m 2 ) 
= J d D l [l 2 - mlnik - hf - m\r[(l - h 



m 



2-1 «3 



(34) 



performing the method presented in Section 1, we obtain the sixty nine non-trivial inde- 
pendent solutions, all analytically continued. We select the convenient solution (consider 
03 = «1 + "2 + «3 + D/2 ) 



= J(3)(£) 

, ot x , a 2 , a 3 , h, k, m , mi, m 2 ) 
= 7r D / 2 (-T^) os (D/2) <n+tta (-Q t3 )_ oa x ^ 3 



-0-3, -ai, — a 2 

1 - (7 3 + « 3 , D/2 



Z| l 2 1 



12 



7715 m 



2 ' 2 ' 2 ' 2 ' 2 
m 2 m 2 m 2 m 2 m 2 



(35) 



where ^3 is a hypergeometric function, that can be given in terms of the generalized Lau- 
ricella functions P7j, (see appendix A). With this solution we can obtain some cases of 
physical interest. The case m 7^ and m\ = m 2 = m, that is applied for example in Higgs 
— ► 77 decay, is given by 

,a 1 ,a 2 ,a 3 ,l 1 ,l 2 ,m ,m) 



_cr 3)ai+D/2 



xR 
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-cr 3 , -«i, -a 3 , -a> 2 



2 > 



2 

12_ 

9 ' 9 



72 
'2 



(36) 



D/2, -ot 2 - a 3 , 1 - ai - D/2 
where i?n is hypergeometric-type function given in the Appendix A. If we make m = and 
mi 7^ m 2 7^ 0, in (p5|), we have the integral that used also in the H decay. The expression 
to this solution is 

j(3) = j( 3 )(D,a 1 ,a 2 ,a 3} li } l 2} m 1 ,m 2 ) 
= n D ^(-mlr(D/2) aia2 (-a 3 U 2 



xR 6 



-cr 3 , -a 2 , -«i 
D/2, 1 - a 3 + a 3 



If m\ If 2 _ l 2 

2' 2' 2' 2 
m 2 777-2 m 2 m 2 



(37) 
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where also Rq is a hypergeometric function, that can be given by generalized Lauricella 
functions, expressed in the Appendix A. The others cases to one-loop three-point function, 
namely: 1) m = and mi = m 2 ; 2) m = mi = m 2 7^ 0; 3)m = m and mi = m 2 = m; 4) 
m = m and m\ = m 2 = 0; 5) m = m 1 = m 2 = 0, are studied in KL . 



VI. FOUR-POINT FUNCTION 

Usually four-point integrals are the most complicated in quantum field theory courses. 
They represent the scattering (2 — > 2) and for this reason are very important in 
phenomenology p0| . In a previous work[[F] two of us studied such integrals — the ones 
that contribute to photon-photon scattering in QED — and did show several hypergeomet- 
ric series representing it. Two of them were calculated before, using Mellin-Barnes approach, 
by Davydychev ||3 1|| , the functions given by Appel's F3 and T,F 2 , the first one a single hyper- 
geometric function and the second a sum of four ones. 

The integral associated to one-loop four-point function, case n = 4, is given 

, Oil) Oi4i lli ^3? m 0t m lt "^2? 



d%[l 2 - mnio - hf - mln(l - l 2 f - mni - h) 



m 



3lQ4 



(38) 



and is far more general than that we considered in [0 and Duplancic and Nizic presented 
in [ 30 1 . This integral can be evaluated by method described in the Section 1 and we obtain 



one thousand and twelve non trivial solution all analytically continued, as can be read 
from table-0. We choose again only the convenient solution, it is written as (consider 04 = 
cti + a 2 + a 3 + a 4 + D/2) 

j(4) = j^^D.a^a^a^a^l^l^l^mQ^i^^m^) 
D/2( -a±)- ai+ai (D/2) a ^ ai _ D/2 (-m^ 

-(j 4 , —aii, -«2, —0C3 
D/2, 1 - 04 + a 4 

(39) 



7T 



-l 2 



4 2 . 



'3 



/2 



'13 



V 2Z mim\ 



ml 



m 



2 ' 



2 ' 2 ' 
mi mi 



m 



2 ' 2 ' 2 ' 2 ' 2 
3 m| 777.3 m 3 m 3 m 3 



2 ' 



where $ is a hypergeometric function (see Appendix A), that can be written in terms of 
generalized Lauricella functions. We observe that above hypergeometric series does not allow 
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one to take limit of vanishing since this it is not defined on that kinematic region. In 
the above solution is possible take out the limits: mo = mi = m 2 = 0; mo = vii\ = and 
m 2 7^ 0; m = and rri\ = m 2 or m 1 7^ m 2 7^ 0; m = mi = m 2 = m 3 = m; and on-shell 
cases, see also J7L [51]. To the special case where m = mi = m 2 = m 3 = we obtain the 
solution 



J(4) = j(4)p 

, Qfi, a 2 , CK3, «4, ^1, ^2, ^3; 



7T 



-Oil, 



-a 2 



(^4)z>/2 

C"4, -"2, -«3 

1 — (T4 + «l , 1 — CT4 + «4 



2 ;2 72 /2 72 

1 2 ^ ^ 



72 72 

L 12 



I 2 



_ _13. ^23 
' 72 ' 72 ' 72 ' 72 

t"3 to to 60 



(40) 



VII. GENERATING FUNCTIONAL FOR PENTAGON INTEGRAL 



Pentagon integrals were studied by Melrose 1 3^] and recently by Bern [[33] and Weinzierl 
and Kosower||34|| , and several authors f35 |, in order to study scattering process where 2 
particles go to 3 particles. 

Let the external legs of the pentagon be (pi, P2—P1, P3—P2, P4 — P3, Pi ), see figure, and 
let us consider firstly the on-shell case. The generating functional for the scalar negative- 
dimensional integrals is written as a special case of the general functional of section 1, 



G° N = I d D q exp ~aq 2 - (3(q + pi) 2 - 7(9 + p 2 ) 2 - % + p 3 ) 2 - co(q + Pi ) 2 



(41) 

then completing the square one can easily integrate and simplify several terms. Eventually 
one get, 



qON 



D/2 



exp 



a^p\ + a6pl + P6sis + (3(f)Su + i4>s 2 a 



(42) 



where we define A = a + (3 + 7 + 8 + u, and three of the Mandelstam variables for the on-shell 
pentagon, 

S13 = (Pi - P3) 2 , s u = (Pi - Pi) 2 , s 24 = (P2 - Pi) 2 - (43) 

Following the usual procedure applied in NDIM we count the total number of solutions we 
will have to deal with. Exponential gives us five sums (Taylor expansion for each argument), 
A other five (multinomial expansion), and the equations are the number of propagators plus 



11 



one, six. Then, there are Cio,6 = 210 possible ways to solve the 6x6 system. The ones that 
have null determinant do not have solution, they are 85, the remaining 125 after properly 
solved will provide us hypergeometric series representing the original Feynman integral. 
On the other hand, projecting out powers of the exponential in fl4"I|), 

00 ( —\ ^01+02+03+04+05 ai aoti a^na^ a 5 

^5=2^ i — i — j — j — j 3* (01,02, 03,04,05), (44) 

Qi ^ 5=0 ai!a 2 !a3!o4!o5! 

where 

J 5 OJV (a 1 ,a 2 ,a 3 ,a4,a 5 )= / A (q 2 ) a ^q + Pl ) 2a *(q + P2? a3 (q + Ps) 2a4 (q + P^ 5 , (45) 



is the negative- dimensional integral. 

We will present two kinds of hypergeometric functions since several others can be written 
immediately if one observe the symmetries of the diagram (or alternatively the symmetries 
of the generating functional), 



{« 3 <-> a 4 ,pi <-> p 4 ~P3,P2 ~Pl ^ P4} ,{ a 5 <-> O4, Pi ^£>3 ~P2,P4 <~> £>2 ~Pl}, 

{ai ^ a 6 , Pi <-> p 3 ~ P2,P2 ~ Pi ^ Pi ~ P3} , (46) 

and so forth, there are several other symmetries of the integral fl4"2p. When we have a hyper- 
geometric series one can obtain some others merely interchanging exponents of propagators 
and external momenta, where three such sets are given in the above list. 



A. Hypergeometric functions for massless on-shell pentagon 

Let us define 05 = 012345 +D/2 and the shorthand notation we will use hereafter in order 
to have a more compact notation, 



Jab Ja T Jb 

and so on. Observe that we use this compact notation only for sum index. Mandelstam 
variables are given by fl43f) . 

The first hypergeometric function representing the Feynman loop integral is given 

-oti, —a 2 , -a 5 , -0-5 - D/2 
l + a 3 -a 5 ,l + a 5 -a 5 



Js ON ({ai}) = n^ ( 4 



(1) c(5) 



P| __pf_ £14 _£13 

S24 S24 S24 S24 



(47) 
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the possible poles are contained in the factor 



if) = 7T D / 2 ( 



-a 3 \a 5 ){-a 5 \a 5 ){cr 5 + D/2\ - 2<r 5 - D/2)s 



as 

21: 



(48) 



where the Pochhammer symbol is given by eq.(21). In the appendix we define all hyperge- 
ometric functions we use in this paper. 

The above hypergeometric function is the "symmetric" solution calculated by 
Davydychev0, using Mellin-Barnes integral representation approach. 

However, in the negative-dimensional approach we obtain, in general, several hyperge- 
ometric functions, which represent it and are related by analytic continuation (direct or 
indirect). We write down other hypergeometric function for the massless scalar pentagon, 

-ax, -a 3 , -a 4 , D/2 + a 123 
1 + a 123 + D/2, 1 - a 145 - D/2 



J 5 ° N ({^}) 



r (2) T -(5) 
1 5 2 4 



1 9 
P2 PI 



£14 £14 
$24 ^13 5 24 ^13 



(49) 



where 



.(2) 



<K D / 2 (-a 2 \a b - k){-a 5 \a 5 - m)(a 5 + D/2\ - 2a 5 - D/2 + a u ) 



cr 5 -a 3 4 
b 14 

-Qt4 e -03 ' 
A 13 6 24 



(50) 



the same set of permutations of exponents of propagators and external momenta can be 
applied to the above result, in order to generate several others. 



B. Hypergeometric functions for massive on-shell pentagon 

We now turn to the case where the five propagators have masses, and let them to be 
distinct. We will present a hypergeometric function that allows us to consider interesting 
particular cases: 2, 3, 4 and 5 equal masses. 

When we deal with massive propagators the system we must solve is greater than the 
related to the massless diagram. Let the masses be m\ attached to the propagators labelled 
as «i, 777J attached to the one labelled as a 2 and so on. 

The generating functional is simply the massless one ( fl~2]) times the exponential containing 
the masses, 

Ggf = G° N exp (am\ + (3m 2 2 + + Qm\ + urn 2 ) , (51) 

the total number of systems we must solve now is given by the combinatorics C15.6 = 5005. 
Among them we pick the most convenient one, i.e., the hypergeometric series that will allow 
us to study several important special cases. We choose to present the result 
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J 5 ON (rn>) 



r (mass)/ ^ 

1 5 (a5,m 5 j^. 



9 5) 



-a 1: —0C 2 , -«3, -OCi, —&5 

D/2, 1 + a 5 - <x 5 



(52) 



where 



^8 



_ Pi 


Z2 = 


A 


2 ' 
















_ 






2 ' 


Z9 = 




mi 




ml 



2.3 



$14 _ £13 

2' ^4 2' ^5 



^24 



m 



2 • 



26 



??? 



2 • 



27 



2 • 



(53) 



which is valid for five different masses, such that ^ and we define 

-,(mass) . m ^ _ D/2 



rr>5; m 5 ) = 7r^(-a 5 \c7 5 )(D/2\ -a,- D/2) (mi) 



2\<T 6 



(54) 



observe that the above 9-fold hypergeometric series allows us to study the cases where any 
of the masses mi,rri2, fn 3 , m^ vanish, or any of the cases where they are equal to m$ (which 
in this kinematical region can not vanish). 

We proceed to special cases now. First, let mi = m 5 . Then, the hypergeometric series in 
the l\ index can be rewritten as (-- - 1 1) - Using Gauss' summation formula |27fl, 

T(c)T(c-a-b) 



2 Fi(o,6,c|l) , 
T(c — a)T(c — b) 

we can exactly sum up the series. Rewriting the Pochhammer symbols involving 1% as, 

X ' "Oil I J12 + — 0" 5 |j'i2345 + ^1234) ( — Oi I J12) ( — 0-5 |ji2345 + few) 



(55) 



ii=0 



/i!(l + a 5 - cr 5 |ji24 + £1234) 



*£ 



1 + a 5 - o<o\ji2i + Z234) 

o"5 + J12345 + ^234^i)(-ai + jn\h) 



(56) 



00 



Zi!(l + a 5 - CT 5 + jl24 + hu\k) 

where we used the property of Pochhammer symbols, 

(a\b + c) = (a\b)(a + b\c), 
thus the series in li is recast as a 2 Fi of unity argument, and then can be summed, 

£ 



(57) 



where 



(I)-" 


(-aib'uX- 


J12345 + ^234) (^4i 1 


-ji235)r(i + a 5 -^)r(A 1 ) 


(I)-. 


(A 2 


- .735) (^3 |j4 + ^234 ) 


T(A 2 )T(A 3 ) 




At = 1 + «i 5 , 


A 2 = 1 + «5, 


A 3 = 1 + a 15 - 0-5, 
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(58) 



(59) 



the same procedure of summing up series of hypergeometric type was applied in [[36 . 

Substituting the above result in fl5"2] ) one obtain the result for the pentagon integral 
having four different masses is given by, 



T ON (m 2 -A\ — r( maSS V™ -m 



-a 2 , -a 3 , -a 4 , -a 5 , -a 5 
D/2, -«i5, 1 + «i 5 - cr 5 



Zi; ...; z 8 



(60) 



where 



y\ v\ _£u _ £13 _£24 

*1 2 ' ^2 2 ' ^3 2 ' ^4 2 ' ^5 2 ' 



mi 



m 
rn 



mi 



m" A 



m 



2 • 



(61) 



where the argument "4" is to remember one the number of different masses and the factor 
involving gamma functions, given by (|54D , is greatly simplified, and we perform the analytic 
continuation in the final step of the calculation (sum the series and then analytically continue 
the Pochhammer symbols). 

The same algebraic manipulations can be made to obtain the cases where the pentagon 
has three different masses ms,m 4 , 7723, 



jfV?;3) = r( mass) (a 12 5;m 5 )^ 



(5) 



-ai, -a 2 , —a 3 , -«4, -a 5 , — cr 5 

D/2, -QJ125, 1 + «125 - 0"5 



where the variables are, 



vl 

z l — o' Z 2 

mi 



vl _ «14 

2 ' Z 3 — 2 ' ^ 4 

mi, mi, 



Sl3 



•§24 



rn 



Zi, z 7 



(62) 



m% m% 
When only two masses are different we can sum the 1% index, 



o> Z § ~~ 21 Zl ~ 2 - (^) 



rn 



mi 



J° N (ml,2) = T { r s \a 1235 ;m 5 )S { 6 5) 



— ai, — a 2 , -a 3 , — a 4 , — a 5 , -cr 5 
D/2, -01235, 1 + «1235 - a 5 



Zi, Zq 



where its six variables are, 



_ vl vl 

Z\ — 2 5 ^2 — o' Z 3 

mi, mi, 



S14 



S13 



■S24 



3" > *4 o ) ^5 9 ) Z 6 9 ) 

mi, mi, mi, 



(64) 



(65) 



finally, the very special case where all five masses are equal, 



J 5 ° N (ml, 1) 



7r ' \m 



-a 5 + D/2\-D/2) 



x<S, 



(5) 



-Oil, -«2, -«3, -«4, -«5, -0"5 
— Q! 12345 



2!i; ...; Z5 



(66) 
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where, 

V\ _ V% S14 _ S13 S 2 4 (af? s 

Z\ — 25 Z 2 — o' Z 3 — 2' Z/l ~~ 2' Z 5 — o' l D ' 

mi, mi, mi, mi, mi, 

observe that two Pochhammer symbols cancelled. The above result is a generalization of 
our previous study of box integrals pertaining to photon-photon scatering0. 



C. Hypergeometric functions for massless off-shell pentagon 



Massless pentagon integrals were necessary to Bern, Dixon and Kosower [j33|l in order to 
study 2 — > 3 scattering, such as e + e~ — > 3 jets. In that work two of the external momenta 
were considered to be massive and the remaining ones massless. 

Here we will present the whole five possibilities: from one to five massive external parti- 
cles. In our approach the difficulty to carry out the integrals with different masses or equal 
ones is the same. So, we prefer the most general case, five distinct masses. 

Hereafter we consider massive external particles, 



V \ = M\ ( Pl -p 2 ) 2 = M 2 l2 , (p 2 -p 3 ) 2 = Ml 3 , ( P3 -p A ) 2 = M 2 A , p\ = Ml (68) 



1. One external leg off-shell 

We present two four-fold hypergeometric series representing the scalar pentagon integral 
with one massive external particle. 



({«*}) = r 5 ( ai ,a 2 ;M 2 )^ 5 
where its five variables are, 



-a 3 , — qj 4 , — a 5 , -ct 5 
1 + a\ — (75, 1 + a 2 — a 5 



Vl V% _ j>13_ _ £14 _ S24 



where the superscript (1) means one leg off-shell and 



(69) 



(70) 



r 5 (ai, a 2 ; M 2 ) = 7 r z? / 2 (-a 1 |a 5 )(-a 2 |a 5 )(a 5 + D/2\ - 2a 5 - D/2){M 



(71) 



observe that the above hypergeometric series does not allow one to take the limit of vanishing 
M since this it is not defined on that point. 
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However, there are others multiple series that allow us to take such limit. One example 

is. 



(5-W/) 



— ai, -a 3 , -a 5 , -cr 5 

1 + a 2 - cr 5 , 1 + a 4 - a 5 



zi; ...; z$ 



define, 



p\ P% S14 s 24 _ M 2 

2l — j ^2 — j Z3 — , Z4 — , Z5 — 

Sl3 «13 «13 «13 «13 



(72) 



(73) 



in the limit of M = the series on j 6 index reduces to its first term, unity, so we are left 
with a four-fold hypergeometric series. 

2. Two external legs off-shell 

Next we put another external leg off mass shell, namely, one consider (pi — P2) 2 = Mf 2 - 
The first hypergeometric series we collect from the total 924 ones (see table [TT]) , is 



J75 (2) (ai,a2,a3,a4,a5) = r 5 (oti, a 2 \ M 2 )Sf 2o//) 



its six variables are given by, 



-CK3, -a 4 , -a s , -0-5 
1 + a\ - cr 5 , 1 + a 2 - 0-5 



(74) 



P2 P3 _ S 13 

_ M 2 ' "' 2 JF- T7 2 ' Z ' 



S24 



M 2 2 



M 2 ' 25 M 2 ' :,i M 2 



(75) 



where the superscript (2) means two legs off-shell. Observe again that the above hyperge- 
ometric series does not allow one to take the limit of vanishing M but admits the limit of 
M 12 = to be taken. 

We can consider also the case where the two masses are equal. Then the hypergeometric 
function in the j-j index reduces to a gaussian one, namely 2-F1, that can be exactly summed 
when its argument is unity. 

Another 6-fold hypergeometric series 



J5 (2) (ai,a 2 , a 3 ,a 4 ,a 5 ) = T' 5 (ai, a 2 , a 3 , ; M 2 2 /pl) 



(5-2o//) 



«2 - o- 5 , -a 2 , -oca, -a 5 , D/2 + a 345 
1 + a 3 - a 5 



Z\\ ...; zq 



(76) 
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where 

_ Vl S13 Si4 S 24 _ M 2 p 2 

Zl - Z2 - M| Z3 'ME' Z4 ~ ~W Zb ~M*2 Z% ~ ~M&' (?7) 
where the factor is defined by, 



2 \ «2 



I* 5 ( ai ,a 2 ,a 3 ,;M? 2 /^) = (- ai \a 5 -a 2 )(-a 3 \a 5 )(a+D/2\-2a 5 -D/2+a 2 ) (^J ( P j)^, 

P2 (78) 

taking M\ 2 = M we can sum up, using Gauss' summation formula[|7]], the j§ series and get, 



Jf\a>i, a 2 ,a 3 ,a A ,a 5 ) = T M T' 5 (ai, a 2 , a> 3 , ; M 2 /p: 



2 /J2\ 



where 



and 



_ vl 

P2 M 1 



(5-2o//) 



"2-^5, -"2, -«4, — Qt5, -D/2 + «345 
1 + «23 - 0"5, 1 - «45 - 0"5 - £/2 



2i; z 5 



Sl3 «14 
' 2 3 — TT^-? ^4 



M 2 



S24 

pi 



M 2 



M 



-a 3 + cr 5 | - o> 2 )(a b + a 4 5 + -D/2| - a 2 ) 



(79) 



(80) 



(81) 



3. Three external legs off-shell 

Continue to study the off-shell pentagon, now where three external particles are massive, 
the third one being (jp 2 ~ P3) 2 = M 23 , an d one such hypergeometric series representation, 



J5 3 \a>i,a 2 ,a 3} a>4,a 5 ) = T 5 (a 3} a 4 ;M 2 



(5-3o//) 
7 



— ai, — a 2 , -a 5 , -a 5 
1 + a 3 - a 5 , 1 - a 4 - a 5 



zi; ...;3 7 



(82) 



where 



M| 3 ' 



22 



M 2 y 



23 



■S13 

M| 3 ' 



z 4 



■S14 
M 2 3 ' 



^5 



■S24 

M| 3 ' 



M 2 M 2 2 

— 2", ^7 " 



Mi 



23 



(83) 



that also admits one to take two masses to vanish. Observe that M 23 can never be zero. 
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4- Four external legs off-shell 



The last step before we put all five legs off mass shell is to make (p^ — Pa) 2 = M| 4 . Then 
we present two sample series: the first is suitable for equal masses limit but non-vanishing, 



Jt\{^}) = T 5 (a u a 2 ;M 2 )St 4 ° ff) 
where its variables are, 



—a 3 , — Qf4, -"5, —&5 

1 + «i - 05, 1 - a 2 - 05 



vl 



vl 



Sl3 



S14 



M 2 ' Z2 M 2 ' M 2 ' 24 M 2 ' 25 



•524 

M 2 ' 



26 



M 2 ' 



27 



(84) 



Mf 4 



M 2 



M 2 



(85) 



and the second one suitable for vanishing masses limit, 



Jt\{ a i}) = r 5 (a 3 , a 5 ;s u )T 8 



(5-4o//) 



— ai, -a 2 , -«4, —0-5 
1 + a 3 - 05, 1 - ct 5 - 0-5 



being 



_ Vl V% «13 

^1 — , ^ 2 — , Z 3 — , 

■524 $24 ^24 



Z A 



S14 
^24 



M 2 
^24 



26 



*24 



^7 



•S24 



28 



Ml 
S24 



(87) 



or when the masses are such that S24 >> M 2 . 

Note that the structure of the above series is similar but not equal. In the first one, only 
one index (j 9 ) appears five times (in other words, the series in j 9 can be rewritten as a 3F2 
function) the other eight indices appear only three times. The second 9-fold hypergeometric 
series has the indices j'3 and appearing five times, what turns it to be more complicated 
than the former. They (equations (84) and (^) ) are very similar but are not the same 
9-fold series. 



5. Five external legs off-shell 

Now we turn to the general case where all five legs are off-shell. Despite there are a great 
number of hypergeometric 9- fold series, see table [H], we present here just two of them, 



J- 5 (5) ({a 4 }) = r 5 ( ai ,a 5 ;M 4 2 )^ 



-a 2 , -a 3 , -a 4 , -05 
1 + a x - 05, 1 + a 5 - 05 



Zi; z g 
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FIG. 1: Scalar massless pentagon. Labels (i,j,k,l,m) represent exponents of propagators and the 
arrows show momentum flow. In the case where the external legs are off-shell we define p\ = M 2 , 
si2 = M 2 2 , s 23 = M| 3 , s 34 = M| 4 , and p\ = M\. 



its nine variables are given as follows, 



Zl 



z 7 



Z2 



Mr Zs 



A. 7 _in 7 

Mf = M 2 ' 3 
M| 4 



Sl4 

M|' 



z 4 



S24 
M 4 2 ' 



^5 



M 2 



^6 



M| : 



Ml' 



i 4 ivl 4 

other hypergeometric 9-fold series that have two masses in the denominator 

J5 (5) («i ? «2,tt3,«4,tt5) = r' 5 (a 4 ,ai,a 5 ;M 4 /M| 4 ) 



x7 -(5-5o//) 



1 + «i4 - cr 5 , 1-0:5 + 0-5 



Zi; 



where we define, 

Zl = 

z 7 = 



PL 
M 4 2 ' 

M 2 2 



Z2 



Z& 



M 2 ' 



z 3 



zg 



gig 

M| 4 ' 
" Ml' 



z 4 



£14 

M| 4 ' 



^5 



■§24 

M| 4 ' 



^6 



M 2 
Mf 



(89) 



(90) 



(91) 



M| 4 ' 8 M| 4 ' 

The last we present is a hypergeometric series that allows one to take the limit of vanishing 
masses (any of them), 



Js I ({«*}) = r 5(«5,ttl,tt3;P2/ S 24) 
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(5-5o//) 



«! — cr 5 , —ax, —0>2, «4 
1 + a 3 - cr 5 , 1 - ax5 + o"5 



Zi; zq 



(92) 



where 



%1 — ~n, Z2 

vi 



£13 

S24' 



23 



£14 

S24' 



M 2 
P2 



S24 



*6 



^24 



2? 



S24 ' 



M 2 s 2 

JW 4 5 24 

^8 — h~i Z 9 ~~ 2~' 

P2 P2 



(93) 



In the negative- dimensional approach we can select the kinematical region one would like 
to study and then work with hypergeometric series defined on that region. All hypergeometric 
series presented in this paper are exact, there are no approximations. They converge very 



fast as we have verified in []37[], a precision of 20 digits can be achieved very quickly 



VIII. HEXAGON FEYNMAN LOOP INTEGRAL 



Recently Bern, Dixon and Kosower|38H studied amplitudes of process involving six parti- 
cles, namely, e + e~ — > 4partons. Also, QCD corrections to e + e~ — > 4jets were calculated by 
Weinzierl and Kosower[^4|; Binoth, Guillet, Heinrich and Schubert |39| show how to reduce 
the hexagon diagrams in order to study Yukawa model at one-loop level. Their works moti- 
vates us to perform such integrals in general cases, i.e., arbitrary exponents of propagators, 
massive external legs and propagators. 

The general formula we calculated in section I gives us the generating functional, 



G 



Hex 



/ \ D / 2 

(?) 

+9(js 35 )} 



cxp 



i + 



1SX3 + (3<j)Sx4 + f3uS 51 + 70S 24 + 7<^25 

(94) 



define o~q = 0123456 + D/2 and ( = A + <p. Total number of systems Cx5j = 6435, such that, 
2790 of them have solutions, the remaining 3645 do not have interest at all. 
The first series we write is of a kind Davydychev|l], |2j called "symmetric", 



Je \{ a i}) = r 6 (a 4 ,a 6 ;s3 5 )5. 



(6— on) 



-ax, -a 2 , -a 3 , a 5 , -cr 6 

1 + «4 — 0"6, 1 + «6 — °"6 



Zl] 



z 8 



(95) 



21 




FIG. 2: Scalar massless hexagon. Labels (ai, 02, 03, 04, 05, a^) represent exponents of propagators 
and the arrows show momentum flow. In the case where the external legs are off-shell we define 
p\ = M 2 , s 12 = S23 = M| 3 , S34 = M| 4 , s 45 = M| 5 and p| = Mf . 



where the function has eight variables, 



T>1 Pt V\ _ -513 «14 _ «51 S24 

Z\ — , Z2 — , Z3 — , Z4 — , Z$ — , Zq — , Z7 — , 

S35 s 35 S35 S35 S35 S35 S35 



£35 
s 35 



(96) 



the superscript (0) means that no external leg is off mass shell. In the appendix we define all 
hypergeometric functions we use in this paper in terms of generalized Lauricella functions. 
The second one pertain to another kind, another kinematical region, 



^6 (0) ({«i}) = r 6 (a 6 ,a 5 ,a 3 ;s24/s25)^8 {6 ° n) 
define, 



"5 - ^6, —Oil, -012, "4, ~«5 
1 + CV 3 — <7 6 , 1 + Ct 56 - CT 6 



Zi\ zg 



,(97) 



V\ V% Vl S13 _ S14 _ S 5 1 _ S35 

Z\ — , Z2 — , Z3 — , Z\ — , 25 — , Zq — , Z7 — , 

S25 S 2 5 S24 S 2 5 S24 S 2 5 S 2 5 

zs = % (98) 
S24 

where in both hypergeometric series the factor Tq means that in equation (ffl]) we change 
cr 5 by a 6 . 
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A. One massive external leg 



Considering one of the external legs of the hexagon to be off mass shell means that one 
must to deal with almost two times more systems of algebraic equations, see table ([TTI|). 

However, some hypergeometric series follow some pattern and we think it can be gener- 
alized for scalar A-point functions. 

The first one we write does not allow one to take the mass to vanish, 



J^({ ai }) = r 6 ( ai ,a 2 ;Mt)St loff) 



-a 3 , -a 4 , -as, «6, -0"6 
1 + a% — er 6 , 1 + a 2 — cr 6 



Zi] Zg 



(99) 



where 



pi 



p\ 



p\ 



Sl3 



ur Z2 Mr Zs Mr Z4 Mr Z5 w -' '" 



gl4 



^7 



•524 



Aff 



2 ; 



•525 



the superscript (1) means that one external is massive. 

Another one, similar to the second we show for the on-shell case, 



s 35 



(100) 



where it has as the previous one nine variables, 



-a 3 , -a 4 , -a 5 , a 6 , a 6 - a 6 
1 + a 2 — cr 6 , 1 + «i 6 — a 6 



Zi; Zg 



aoi) 



P2 ^3 P\ _ S 13 "514 

Zi_ Mf' Z2 ~wr z "~w ZA ~w Z5 ~w Z6 

_ S35 _ Ml 

z 8 — j Zg — , 

S51 S51 



■524 _ ^25 



(102) 



where r' 6 (.) = r' 5 (.) substituting <r 5 by erg. 



B. Six massive external legs 

Instead of presenting the particular cases - as we have done for the pentagon - where two, 
three, four and five external are massive, we jump toward the most general case, namely, 
six massive external legs, 



J?\{<Xi}) = T 6 (a 1 ,a 6 ;M'£)S : 



(6-60//) 

14 



-a 2 , -« 3 , -"4, -«5, -cr 6 
1 + «i — cr 6 , 1 + a e — (Jq 



Z\\ 



Z\A 



(103) 
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where 





pi 




pi 




Z\ = 


Mr 


Z2 = 


Mr 


z 3 








s 35 




z 8 = 


w 


Z9 = 




ZlO 


Z14 = 


Ml 








Ml 









2 

P4 _ s 13 _ s 14 _ S51 _ s 24 

M? _ _ Mf 2 _ M| 3 _ M| 4 



Mf Zu Mi' Zl2 Mi' Zl3 Mi' 



(104) 



where £7 = J123456789 and Era = 7112345- From the above result we can infer several special 
cases, namely, on-shell external legs, equal masses external legs and so on. To work out 
these particular cases one can proceed on the same way we did in the previous sections and 
subsections. 



C. On-shell hexagon with 6 massive propagators 

In this subsection, the last one, we present for completeness a result for the hexagon 
scalar integral, where all its external legs are massless and on-shell and its propagators have 
distinct masses, Mi,...,M 6 . 

We pick a sample hypergeometric series, that allows us to obtain several particular cases 
(equal masses, or some of them null). The one that is not contained in this kinematical 
region is the case that has M§ = 0. 

Let us call such 14-fold series J§ m 



rr f( 1 \ -r^imass) 1 „, \ 0(6— Gmass) 



-ati, -a 2 , — ot 3 , -a 4 , -a 5 , -cr 6 
1 + a 6 - a*, D/2 



Zi\ Zu 



(105) 



where 

zi 
z 8 = 



Pi 



Z2 



pi 



Mi' 

S35 S35 



Mf 



Z3 



pi 



Mi 



2 • 



Mr 



Z9 



Mi' 



ZlO 



M[ 

w 



zu 





Sl3 




S14 


Za = 


Mi 


, z b = 


Ml' 


Ml 




Ml 




~mI' 


Z\2 


~mI' 


Zl3 = 



S51 



S24 



Ze Ml' Z7 ~~ Mi 



Mi 



M. 



2 ■ 



ZlA 



Ml 
Mi' 



(106) 



where rg mass ^ is given by equation (|54j) changing cr 5 by <j 6 . So we finish our study on six-point 
functions. Special cases where two or more masses are equal can be obtained summing up 
the series, as we did in previous sections. 
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IX. CONCLUSIONS 



Higher one- loop n-point integrals are becoming more important to scattering process |33 
39. 40] in the standard model. We presented in this work a general formula for such integrals, 
and did show how to apply it to several cases of interest, from one through six-point integrals. 
Being its external legs on or off mass shell, its propagators massless or massive, and their 
exponents arbitrary. Working with different masses or equal ones raises the same difficult in 
the NDIM approach, as well as exponents of propagators, taking the most general case does 
not imply in additional technical difficulties. We choose sample hypergeometric solutions, 
the ones in which one can take interesting particular cases. Of course, it is not possible to 
present all of them (see tables I, II and III). The same procedure can be followed in order to 
calculate even higher (N > 7) scalar integrals. 
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APPENDIX A: SOME HYPERGEOMETRIC FUNCTION USED 

We present here some of hypergeometric or hypergeometric-type functions used in this 



iA:B L 



paper all in terms of the generalized Lauricella functions ]41]], that is expressed by 

[a:aW,.-,« W ];[6 (1) :/?«];...; [&W : pW] 
[c:^),..^^];]^:^];...;]^):^)] 



Z\\ zm 



E 

jl,...,jjv=0 



(N) . 



ll i= i(Ci) (D il+ .„ + <*> iw 



D(*0, .(A^)^ 



5 i ON 



Z 



z 



JN 

N 



where 



J = (oi : a\\...,a\ '),...,{a A : oc\',...,a\ ') 



S N ) 



(1) 



(k) 



}(*) 



c : 7 (1) , 7 W ] = (ci : 7*\ tD. («c : 7^ 7^) 



,(!) 



(Al) 
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with k — 1, N. The parameters a, (3, 7, 5 can be non-negative integers in hypergeometric 
functions and too negative integers in hypergeometric-type functions. From (|A1|) we can 
extract all the functions used in this paper, that is, 



Re, — Re 



Xi 


X2,X 3 








z \\ z 2\ Z3; Z A 


£4 


X 5 





-^1;0;0;0;0 



(an : 1, 1, 1, 1), (x 2 : 1, 1, 1, 0), (x 3 : 1, 0, 0, 1) 
(x 4 : 1,0, 1,1), (x 5 : 1,1,0,0) 



z 2] Z 3 \ Za 



R 



11 



R 



11 



-^3;0 



Z\ \ ...; Z4 



Xi,X2, X 3 , X4 
X5, Xq, X-j 

{x x : 1), (x 2 : 1, 1, 0, 1), (x 3 : 1, 0, 1, 1), (x 4 : 1, 0, 1, 0) 
(x 5 : 1, 0, 1, 1), (x 6 : 1, 0, 2, 1), (x 7 : 0, 1,-1,0) 



z\ ; . . . ; Z4 



where F 3 4 ; : ° = F 3 : : j : and (x l : 1) = (x 1 : 1, 1, 1, 1), 
^3 = ^3 



F 3:0 
r 2;0 



x 1 ,x 2 ,x 3 








zi; . 




X4, X5 







( Xl : 1, 1, 1, 1, 1), (x 2 : 1, 1, 0, 1, 0), {x 3 : 1, 0, 1, 0, 1) 
(x 4 : 1,0, 0,1,1), (x 5 : 1,1,1,0,0) 



where F™ = F™$$, 



$ = $ 



F 4:0 
r 2;0 



Xl, X2, X 3 , £4 
£5, 

(xi : ai), (x 2 : a 2 ), (23 : a 3 )(x 4 : a 4 ) 
(x 5 : a 5 ), (x 6 : a 6 ) 



Z\\ Zg 



(A2) 



(A3) 



(A4) 



(A5) 



where F^r^ = F 2 4 : ; % . ; / ; % ; °. ; and (x 1 : o x ) = {x x : 1,1,1,1,1,1,1,1,1), (x 2 : a 2 ) 
(x 2 : 1,1,0,1,0,0,0,1,0,0), (x 3 : a - 3) = (x 3 : 1,0,0,1,1,0,0,1,0), (x 4 : a 4 ) 
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(x 4 : 0,1,0,1,0,1,0,0,1), (x 5 : a 5 ) = (x 5 : 1,1,1,1,1,1,0,0,0) and (x 6 : a 6 ) = (x 6 : 
1,1,0,1,0,0,0,0,1). 

For the pentagon integrals we define the hypergeometric series, where "5" and A 
represent pentagon and the number of variables the series has, respectively. 
In equation ( pE7j ) the result was written in terms of, 



c(5) _ <j(6) 



F 4:0 



Xi, X2, X3, X 4 
X5 , Xg 

(xi : 1, 1, 0, 0), (x 2 : 0, 0, 1, 1), (x 3 : 0, 1, 0, 1), (x 4 : 1) 
(x 5 :0,l,l,i),(x 6 :l,l,0,l) 



Z\\ ...; Z4 



(A6) 



Equation (|49"D was written as a series defined by, 

(5) _ ^(5) 



-*4 — -*4 



Z\ \ ...; z 4 



7 4:0 
2:0 



Xi, X 2 , X3, X 4 
X 5 ,X 6 

(X! : 1, 1, 0, 0), (x 2 : 1, 0, 1, 0), (x 3 : 0, 1, 0, 1), (x 4 : -1, 0, 0, 1) 
(x 5 :0,0,l,l),(x 6 :0, 1,-1,0) 



21; Z4 



(A7) 



Five variables series was obtained as a result for the massive pentagon in the special case 
where all five masses are equal, 



5 (5) _ £(5) 



^i;0 



3^1)^2) 3^4) ^5; ^6 
• r 7 



Zi; ...; 2g 



(xi : 11000), (x 2 : 00110), (x 3 : 10001), (x 4 : 01010), (x 5 : 00101), (x 6 : 1) 
(x 7 : 2) 



21; ...; z 5 



(A8) 



where we did not write the commas, i.e., (x 4 : 01010) = (x 4 : 0, 1, 0, 1, 0) and so forth. 
Six variables is one of the representation for the pentagon with two distinct masses, 



c(5) 
<->6 



F 6:0 
r 3:0 



Xl,X2, X3, X 4 , X5, Xq 
X7, Xg, X9 

(xi : 110000), (x 2 : 001100), (x 3 : 100010), (x 4 : 010101), (x 5 : 001010), (x 6 : 1) 
(x 7 : 111110), (x 8 : 212120), (x 9 : -1, 0, -1, 0, -1, 1) 
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The reader can observe that the following result and previous two follow a pattern, seven 
variables, for three different masses, 



o(5) _ o(5) 



-^3;0 



Xi, X2, X 3 , X 4 , X5, Xq 

( Xl : 1100000), (x 2 : 0011000), (x 3 : 1000110), (x 4 : 0101001), (a* : 0010100), (x 6 : 1) 
(x 7 : 1111100), (x 8 : 1121100), (x 9 : 00, -1, 0011) 



Zl 



c(5) . . o(5) 

Do — Oo 



^3:0 



Zi; z 8 



Zi, z 8 



The eight-fold series representing the case where the pentagon has four different masses, 

(ai), (a 2 ), (a 3 ), (a 4 ), (a 5 ), (a 6 ) 
(a 7 ), (a 8 ), (a 9 ) 

(ai), (a 2 ), (a 3 ), (a 4 ), (a 5 ), (a 6 ) 
(x 7 : 11111000), (x 8 : 1111100), (x 9 : 00010111) 

(AH) 

where, a x = (x : : 11000000), a 2 = (x 2 : 00110100), a 3 = (x 3 : 10001010), a 4 = (x 4 : 
01010001), a 5 = (x 5 : 00101000), a 6 = (x 6 : 1). 

Finally the nine-fold series, representing the on-shell pentagon with five distinct masses, 

(bi),(b 2 ),(b 3 ),(b 4 ),(b 5 ) 

(x 7 ), (x 8 ) 

(61), (6a), (6s), (64), (65) 
(x 7 : 111110000), (x 8 : 1101011111) 

(A12) 



<j(5) _ <j(5) 



Zi, ...;zg 



— P5:0 

- -^2:0 



Zi, Zg 



where, 61 = (xi : 110001000), 6 2 = (x 2 : 001100100), 6 3 = (x 3 : 100010010), 6 4 = (x 4 : 

010100001), 65 = (x 5 ■ 1). 

The off-shell pentagon starts with a five variables Lauricella function, in the case where 
one leg is off mass shell, 



,(5-lo//) = £(5-lo//) 



Xi, X 2 , X3, x 4 
X5 , Xg 



zi; Z5 



P 4:0 



(xi : 10001), (x 2 : 01100), (x 3 : 00011), (x 4 : 1) 
(x 5 : 01110),(x 6 : 11001) 



Zi, z 5 



(A13) 
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the following function represents eq.( 



jis-ioff) _ rj-{b-loff) 



F 4:0 



zi;...;zs 



x ll %2i x 3i %4 
X5 5 Xg 

(xi : 11001), (>2 : 10010), (x 3 : 00110), (x 4 : 1) 
(x 5 : 11010), (ar 6 : 10111) 



(A14) 



we presented two results when two external legs are off-shell, the first one was given in terms 
of. 



q (5-2off) _ q (5-2off) 



F 4:0 
r 2;0 



and the second one, 



%l, %2, X3, X 4 
X 5 ,X 6 

{x x : 100011), (x 2 : 011000), (x 3 : 000110), (x 4 : 1) 
(x 5 : 001111), (x 6 : 110010) 



Zi, ...; z 6 



(A15) 



rAh-loff) _ rAh-loff) 

*6 — J 6 



P 5:0 



°1; a 2i a 3i a 4i a 5 
Xg 

(ai), (a 2 ), (a 3 ), (a 4 ), (a 5 ) 
(x 6 : 111010) 



Z\] zq 



Z\\ ...; Zq 



(A16) 



where (a x ) = (xi : 10010,-1), (a 2 ) = (x 2 : 010011), (a 3 ) = (x 3 : 110000), (a 4 ) = (x 4 : 
001100), (05) = {x 5 ■ 000, -1, 11), (a 6 ) = (x 6 : 111010). 
The special case of M — M 12 gives rise to, 

X\, X2, X 3 , X 4 , X5 
X 6 , X7 



^(5-2o//) _ y{b-2off) 



z\\ Z5 



F 5:0 
r 2;0 



(xi : 1001, -1), (x 2 : 01101), (x 3 : 11000), (x 4 : 00110), (x 5 : 000, -1, 1) 
(x 6 : 1000,-l),(x 7 : 1111,-1) 



Zi;...;z$ . 
(A17) 



Three external legs off shell was written as a function, 

X\, x 2 , x 3 , x 4 , 

X5, Xg 



<j(5-3o//) _ <j(5-3o//) 

O7 — Oy 



Zi; z 7 
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— P 4: 

- -^2:0 



(xi : 1100010), (x 2 : 0011011), (x 3 : 00011000), (x 4 : 1) 
(x 5 : 0111010), (x 6 : 1001111) 



(A18) 



Four external legs, the first result presented, 



c (5-4o//) _ c (5-4o//) 



- ^4:0 

~~ *2;0 



z\\...;z% 



0,1,0,2, Ct 3 , Cl 4 

(ai), (a 2 ), (a 3 ), (a 4 ) 

(x 5 : 00111111), (x 6 : 11001011) 



Z\, 



(A19) 



where {a x ) = (x 1 : 10001110), (a 2 ) = (x 2 : 01100011), (a 3 ) = (x 3 : 00011001), (a 4 ) = (x 4 : 
1), and the second is written as the generalized Lauricella function, 



j-(5-4o//) = T (5-4o//) 
8 8 


ai, 


a 2 , a 3 , a 4 








x 5 


Xq 







4:0 



-^2;0 



(ai), (a 2 ), (a 3 ), (a 4 ) 

(x 5 : 00111001), (x e : 11101110) 



Zi; z 8 



(A20) 



where (ai) = (xi : 11001000), (a 2 ) = (x 2 : 00111100), (a 3 ) = (x 3 : 01100011), (a 4 ) = (x 4 : 
!)• 

Finally, we define the hypergeometric functions for the five external legs off-shell pen- 
tagon, 



0(5-50//) c(5-5o//) 

o 9 — o 9 


ai, 


a 2 , a 3 , a 4 


zi; . 


.; 




x 5 


x 6 







- -^2:0 



(ai), (a 2 ), (a 3 ), (a 4 ) 

(x 5 : 001110111), (x 6 : 111001110) 



Z\, Zg 



(A21) 



where (ai) = (xi : 001101100), (a 2 ) = (x 2 : 100010110), (a 3 ) = (x 3 : 011000011), (a 4 ) 
(x 4 : 1), the second one, 



^-(5-5o//) ^( 5 -5o//) 


ai, 


a 2 , a 3 , a 4 




•i ^9 




x 5 


x 6 
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P 4:0 
-^2;0 



(ai), (a 2 ), (a 3 ), (a 4 ) 

(x 5 : 0, -1, 011010, -1), (x 6 : 111001110) 



Zi; 2 9 



(A22) 



where (ai) = (xi : 110001001), (a 2 ) = (x 2 : 001101100), (a 3 ) = (x 3 : 100010110), (a 4 ) 
(x 4 : 00111011, -1), and the third, 



w (5-5o//) = Wg (5- 5 °//) 



- /?5:0 
~~ *2;0 



Zi; z 9 



ai, a 2 , a 3 , cz.4, as 
(ai), (a 2 ), (a 3 ), (a 4 ) 

(x 6 : 111100110), (x 7 : 0100110, -1, -1) 



Z\\ £g 



(A23) 



[a 2 



(x 2 : 100100011), (a 3 ) = (x 3 



where (a x ) = (x x : 01101110,-1), 
011110000), (a 4 ) = (x 4 : 110001100) 

The hexagon integrals are also defined in terms of generalized Lauricella functions, the 
firs result we presented for the on-shell case is given as, 



(6— on) 



s> 



-^2;0 



z±; z& 



(6-on) °1' °2) Q 3 , °4, ^5 
^6,^7 

(ai), (a 2 ), (a 3 ), (a 4 ), (a 5 ) 
(x 6 : 10101111), (x 7 : 11111010) 



(A24) 



where (ai) = {x x : 11100000), (a 2 ) = (x 2 : 00011100), (a 3 ) = (x 3 : 10000011), (a 4 ) = (x 4 : 
00101010), (a 5 ) = (x 5 : 1). The second one, 



(6— on) 



r 8 



F 5:0 
r 2;0 



(6-on) 01,02,03,04,05 
^6,^7 

(ai), (a 2 ), (a 3 ), (a 4 ), (a 5 ) 
(x 6 : 01111110), (x 7 : 1101000, -1) 



Zi,...;z$ 



(A25) 



where (ai) = (xi : 1101011,-1), (a 2 ) = (x 2 : 11100000), (a 3 ) = (x 3 : 00011100), (a 4 ) 
(x 4 : 01010010), (a 5 ) = (x 5 : 00101001). 
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When one of the external particles is massive, or an external leg is off mass shell, the 
integral is written in terms of 9-fold Lauricella function, 



q (6-loff) _ <j(6-lo//) 



^2;0 



Z\\ Zg 



ai, a 2 , a 3 , 04, 0,5 

(ai), (a 2 ), (a 3 ), (a 4 ), (a 5 ) 

(x 6 : 00011111), (x 7 : 111000111) 



Zi; zg 



(A26) 



where (ai) = (x 1 : 100000110), (a 2 ) = (x 2 : 010100001), (a 3 ) = (x 3 : 001010100), (a 4 ) = 
(x 4 : 000001011), (a 5 ) = (x 5 : 1), the next we presented is also for the case of one massive 
external particle, 



T-Ce-lo//) _ ^(6-loff) 
2 9 — 2 9 



F 5:0 
r 2;0 



di, d2, ds, d4, 05 

x 6 ,x 7 

(ai), (a 2 ), (a 3 ), (a 4 ), (a 5 ) 

(x 6 : 111001110), (x 7 : 00011100, -1) 



Z\\...\Zg 



(A27) 



where (<n) = (xi : 100001100), (a 2 ) = (x 2 : 010100010), (a 3 ) = (x 3 : 001011000), (a 4 ) = 
(x 4 : 000000111), (a 5 ) = (x 5 : 11111100, -1). 

The most complicated integral, the hexagon with all external legs massive was given as 
a 14-fold series, 



^(6-60//) _ ^(6-60//) 



5:0 



Z\] ...] Z14 



— f? 0:lJ 

- -^2;0 



ai, a 2 , a 3 , a 4 , 05 
x 6 ,x 7 

(ai), (a 2 ), (a 3 ), (a 4 ), (a 5 ) 

(x 6 : 00011111101111), (x 7 : 11111010011110) 



Z\\...\Z\4 



(A28) 



where (a x ) = (xi : 00011111000), (a 2 ) = (x 2 : 10000011001100), (a 3 ) = (x 3 : 
01010000100110), (a 4 ) = (x 4 : 00101010000011), (a 5 ) = (x 5 : 1). 

Finally, the last generalized Lauricella function, for on-shell hexagon with six different 
massive propagators, 



r.(6— dmass) c(6— Gmass) 

o 14 — o 14 



di, a 2 , a 3 , a 4 , 05, a6 



Zi; ...; Zu 
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— P6:0 

- -^2;0 



(ai), (a 2 ), (a 3 ), (a 4 ), (a 5 ), (a 6 ) 

(x 7 : 11111110011111), (x 8 : 11111111100000) 



z\ \ ...; Z14 



(A29) 



where (ai) = (xi : 11100000010000), (a 2 ) = (x 2 : 00011100001000), (a 3 ) = (x 3 : 
10000011000100), (a 4 ) = (x 4 : 01010000100010), (a 5 ) = (x 5 : 00101010000001), (a 6 ) = 
(a* : 1). 
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TABLE I: Off-shell Scalar Box Integral: number of systems, solutions and kind of series 


Off-shell Box 

V / ±L Oil v.- ± L 1 r Wjl 


Total 


Sol iit,i oris 


TVo solution 


T~rvnprp"POTTi sprips 


Massless 


252 


162 


90 


5-fold 


1 mass 


462 


267 


195 


6-fold 


2 masses 


792 


426 


366 


7-fold 


3 masses 


1287 


663 


624 


8-fold 


4 masses 


2002 


1012 


990 


9-fold 


TABLE II: 


Pentagon 


Integral: number of s) 


r stems, solutions 


and kind of series 


Pprit n crrm 

J. \~>LL L cl ii w±± 


Total 


Slolntiorm 


TVn solution 


l-Tv'nprp'pn'm series 


On-shpll 

N / -LJ. OllV'll 


210 


125 


85 


4-fold 


1 off-shell 


462 


247 


215 


5-fold 


2 off-shell 


924 


474 


450 


6-fold 


3 off-shell 


1716 


855 


861 


7-fold 


4 off-shell 


3003 


1518 


1485 


8-fold 


5 off-shell 


5005 


2530 


2475 


9-fold 


TABLE III: Hexagon 


Integral: number of systems, solutions and kind of series 


Hexagon 


Total 


Solutions 


No solution 


Hypergeom. series 


On-shell 


6435 


2790 


3645 


8-fold 


1 off-shell 


11440 


4736 


6704 


9-fold 


2 off-shell 


19448 


7155 


12293 


10-fold 


3 off-shell 


31824 


12408 


19416 


11-fold 


4 off-shell 


50388 


19484 


30904 


12-fold 


5 off-shell 


77520 


30410 


47110 


13-fold 


6 off-shell 


116280 


45615 


70665 


14-fold 



